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, Abstract. We investigate the gravitational wave background produced by mag- 

netars. The statistical properties of these highly magnetized stars were derived 
by population synthesis methods and assumed to be also representative of extra- 

OO 

galactic objects. The adopted ellipticity was calculated from relativistic models 
using equations of state and assumptions concerning the distribution of currents 
in the neutron star interior. The maximum amplitude occurs around 1.2 kHz, 
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corresponding to Q. gw ~ 10 9 for a type I superconducting neutron star model. 
The expected signal is a continuous background that could mask the cosmological 
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5^ ■ contribution produced in the early stage of the Universe, 
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H 1. Introduction 



The stochastic background is a rich source of gravitational waves (GW). The detection of 
a cosmological gravitational background (CGB) is an unique way to probe the very early 
stages of the universe up to the limits of the Plank era and the big bang (Grishchuk et al., 
2001). In addition to the CGB, the emission from a large number of unresolved sources is 
expected that produces the astrophysical gravitational background (AGB). Because the 
AGB mask the relic signal, the knowledge of its properties such as the spectral energy 
distribution is of fundamental importance to determine the best frequency domain in 
which to search for the CGB in order to optimize the subtraction techniques. A large 
number of astrophysical processes able to generate a stochastic background have been 
investigated. On the one hand, distorted black holes (Ferrari et al. 1999a; de Araiijo et al. 
2000) and bar mode emission from young neutron stars (Regimbau 2001) are examples of 
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sources able to generate a shot noise signal (the time interval between events is large in 
comparison to the duration of a single event), while supernovas (Blair et al. 1997; Coward 
et al. 2000; Buonanno et al. 2004) are expected to produce an intermediate "popcorn" 
noise. On the other hand, the contribution of tri-axial rotating neutron stars results in 
a truly continuous background that may compete with the relic one in the frequency 
domain of ground based interferometers (Regimbau & de Freitas Pacheco 2001b, here- 
after RPOlb). Neutron stars (NS) are one of the most promising sources for observation 
by ground based gravitational wave interferometers such as LIGO, VIRGO, GEO and 
TAMA (see de Freitas Pacheco 2001 for a recent review). Various processes related to NS 
able to produce large numbers of gravitational waves (GWs) have been investigated: free 
precession of an axisymmetric star induced by a misalignment between the spin and the 
symmetry axes (Zimmcrmann & Szedenits 1979), r- or bar-mode instabilities excited by 
the so-called CFS (Chandrasekhar-Fricdman-Schultz) mechanism, which may occur in 
fast and hot rotating neutron stars (Andersson & Kokkotas 2001 and references therein) 
and mini gravitational collapse induced by a phase transition in the core (Marranghcllo, 
Vasconcellos & de Freitas Pacheco 2002). However, the emission produced by a tri-axial 
rotating star is still one of the most studied mechanisms. Different scenarios leading to 
a distorted star have been discussed in the literature, as anisotropic stresses from strong 
magnetic fields and tilting of the symmetry axis during the initial cooling phase when 
the crust solidifies. Bildsten (1998) pointed out that a neutron star in a state of accretion 
may develop non-axisymmetric temperature variations on the surface, which produce 
horizontal density patterns able to create a large mass quadrupole moment, if the elastic 
response of the crust is neglected. More detailed calculations (Ushomirsky, Cutlerand 
Bildsten 2000) indicate that the inclusion of the crustal elasticity decreases the expected 
mass quadrupole by a factor of 20-50, reducing considerably the predicted GW emis- 
sion. More recently, elastic deformations of compact objects with "exotic" equations of 
state have been considered by Owen (2005). He showed that solid strange stars could 
sustain ellipticities as high as few times 10~ 4 , considerably higher than estimated values 
for conventional neutron stars (Thorne 1980). The distortion induced by magnetic fields 
becomes significant in highly magnetized neutron stars, overwhelming the "flattening" 
due to a fast rotation. The existence of neutron stars with magnetic fields in excess of 
10 14 G (called magnetars) was firstly suggested by Thompson & Duncan (1992). They 
have shown that in fast newborn neutron stars, the dynamo mechanism could generate 
magnetic field strengths up to 10 16 G. The existence of magnetars is supported by the 
observation of "soft gamma repeaters" (SGRs) and anomalous X-ray pulsars (AXPs), 
whose rotation periods and deceleration rates, if interpreted in terms of the canonical 
magnetic dipole model, suggest that these objects are associated with young highly mag- 
netized NS (Kouveliotou et al. 1998, 1999; de Freitas Pacheco 1998; Heyl & Kulkarny 
1998; Mereghetti 1999). On the one hand, magnetars may be considered as objects hav- 
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ing an origin and evolutionary path different from "classical" radio-pulsars but, on the 
other hand, simulations based on population synthesis methods suggest that NS are born 
with a large variety of rotational periods, magnetic fields and that magnetars are objects 
simply born in the high end of the magnetic field distribution (Regimbau & de Freitas 
Pacheco 2001a, hereafter RPOla). According to these simulations, a non-negligible frac- 
tion of NS are born with magnetic fields higher than 10 14 G, raising the possibility that 
these objects may also generate a continuous stochastic background. In RPOlb, the con- 
tribution of tri-axial rotating neutron stars to the gravitational wave background was 
calculated under the assumption that the initial period and magnetic field distributions 
derived from population synthesis of galactic pulsars by Regimbau & de Freitas Pacheco 
(2000, hereafter RP00) could be extended to the extragalactic population. In this paper, 
we adopt a similar procedure to compute the contribution of magnetars to the GW back- 
ground, but with a fundamental difference. In RP00 the equatorial ellipticity was taken 
as a free parameter, whereas in the present work the equatorial ellipticity is estimated 
from NS models including different magnetic field configurations and equations of state. 
This paper is organized as follows: in Section 2 we discuss the properties of the magne- 
tar population, including the magnetic field distribution and the equatorial deformation, 
in Section 3 their contribution to the gravitational wave background is calculated, in 
Section 4 the expected signal is compared to the CGBs and the detection possibility by 
present and future generations of GW detectors is discussed and in Section 5 the main 
conclusions are summarized. 

2. Properties of magnetars 

2.1. The magnetic field distribution 

In our previous work (RP00), the population synthesis method was used to recover the 
statistical properties of the real population of radio pulsars, since most of the objects 
are hidden by strong selection effects, such as the flux density limit of radio telescopes 
and their orientation with respect to the pulsar emission beam. These simulations were 
recently upgraded in order to take into account the new pulsars discovered by the Parkes 
multibcam survey In agreement with our previous analysis, the initial period follows 
a normal distribution of mean 240 ms with a dispersion of 80 ms, while the magnetic 
field has log-normal distribution of mean log B = 13.0 (in Gauss) with a dispersion of 
&iogB —0.8. For details of these simulations, the reader is referred to RP00. 

As shown in Fig. 1, the magnetic field distribution of the observed population matches 
quite well that derived from the observed period deceleration rate and the "canonical" 
magnetic dipole model. In this case, the mean field of the observed population is about 
1.6 x 10 12 G and the distribution has an apparent cutoff around B = 10 14 G. As we have 
already stated (RPOla), the average magnetic field of the real (or "hidden") population 
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Fig. 1. Comparison between the simulated and the observed distribution (derived from 
the "canonical" magnetic dipole model) of the magnetic field (in Gauss). 



is one order of magnitude higher and the fraction of magnetars, defined as objects with 
fields above 10 14 G, is around 8%. A comparison between the magnetic field distributions 
of the observed and the "unseen" populations is given in Fig. 2. Notice that objects with 
fields as high as 3 x 10 15 are expected to exist in the Galaxy. 

2.2. Deformation by magnetic fields 

The deformation of magnetized Newtonian stars was already discussed in the early fifties 
by Chandrasekhar & Fermi (1953) and Ferraro (1954). The gravitational wave emission 
from magnetic distorted stars was considered by Gal'tsov & Tsvetkov (1984), Bonazzola 
& Gourgoulhon (1996), Konno, Obata & Kojima (1999, 2000), Palomba (2001), among 
others. The deformation of a slowly rotating magnetized star can be expressed as the sum 
of the contribution of three main terms: the first and the most important corresponds 
to the Lorentz force, induced by current flows in the highly conductive NS interior; the 
second term represents variations of the gravitational potential, a consequence of the 
distortion itself and the third is a purely relativistic term arising from the definition of 
the circumferential radius R = Lj (2tt), where L is the length of the equator as measured 
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Fig. 2. Magnetic field distributions of the observed (dashed line) and the unseen popu- 
lation (full line) of pulsars 

by a non-rotating observer. In general, the ellipticity induced by magnetic field effects 
can be expressed by the dimensionless ratio (Konno, Obata & Kojima 2000) 

B 2 R A 

e b = g-^^sm 2 a = 1.9 x lO^gB^R^M^jsi^a (1) 

where B, R and M are respectively the magnetic field at the surface, the radius and the 
mass of the star, a is the angle between the spin and magnetic dipole axes, while g is 
a dimensionless parameter depending on both the equation of state (EOS) and on the 
magnetic field geometry. It is equivalent to the parameter (3 introduced by Bonazzola & 
Gourgoulhon (1996). 

In the case of an incompressible fluid star with a dipole magnetic field, considered 
by Ferraro (1954), the deformation parameter is g = 12.5. Relativistic models based on 
a polytropic EOS, e.g., P oc p v and a dipole field geometry give similar values (Konno, 
Obata & Kojima 2000). NS models built with the EOS UV14+TNI (Wiringa, Fiks & 
Fabrocini 1988) and with non-superconducting neutron star matter lead to deformation 
parameters g of the order of the unity for currents concentrated in the crust and a few 
times higher if currents are concentrated in the core (Bonazzola & Gourgoulhon 1996). 
These values correspond to ellipticities of about (2 — 20) x 10~ 8 for fields typically of the 
order of 10 14 G, below the limits estimated by RPOlb to produce a detectable signal by 
the present generation of interferometric detectors. A more favourable situation occurs if 
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Fig. 3. Simulated distribution probability for the expected magnetic ellipticity of mag- 
netars 

neutron stars have a superconducting interior of type I. In this case, the magnetic field 
permeates only the very outer layers of the star and distortion parameters as high as 
g w 520 can be obtained (Bonazzola & Gourgoulhon 1996). From eq. 1 and the magnetic 
field distribution of the actual population derived in the previous section, the ellipticity 
distribution of magnetars can be calculated (Fig. 3). 

Our simulated data can be quite well fitted by the function 

p(x) = 2.36 exp(-0.2(x + 9.7) 2 ) if x > -7.7 (2) 

where x = logfe^) — 21og(g) 

According to our simulations, even if magnetars with ellipticities as high as ~ 0.01 
are expected to exist, their gravitational radiation remains small, since most of them 
have low rotation frequencies (RPOla). This implies that in spite of their high ellipticity, 
these objects are always decelerated by the magnetic braking mechanism, a point to be 
retained when the spectral energy is computed (see Section 3). To illustrate this point, 
the distribution of the ratio R between the two energy loss mechanisms, namely 

R = ^MUL = 1.36 x lQ- 12 (^-) 2 B 2 14 sin 2 aP- 2 (3) 

is shown in Fig. 4. The numerical value was calculated for a "canonical" NS of mass M = 
1.4 M , radius R = 10 km, moment of inertia I = 1.4x 10 45 gem 2 and B = 10 14 G. Notice 
that even for "extreme" cases, e.g., magnetars in the tail of the period and magnetic field 
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Fig. 4. The ratio between the gravitational and magnetic energy loss rates for magnetars 

distributions, the aforementioned ratio is R ~ 0.002 << 1 for objects with P <~ 1ms, 
B 14 , ~ 1 and g ~ 520. 

3. The stochastic gravitational wave background 

3.1. Spectral properties 

The superposition of the gravitational wave emission from different objects and in par- 
ticular from magnetars formed since the beginning of the star formation in the universe 
produces a stochastic background whose detection might be possible by employing cross 
correlation techniques between two antennas (Grishchuk 1976; Christensen, 1997). 

The spectral properties of the stochastic background are characterized by the dimen- 
sionless parameter (Ferrari et al. 1999): 

ttgw(Vo) = -^—V F Vo (4) 
C Pc 

where v Q is the wave frequency in the observer frame, p c is the critical mass density 
needed to close the Universe, related to the Hubble parameter H by 

F Vo is the gravitational wave flux at the observer frequency v integrated over all sources 
defined by 

F Vo = f™" UJR(z) (6) 
Jo 
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where f Vo is the gravitational flux of a single source located between z, z + dz and dR(z) 
is the source formation rate. Both quantities are given by 

f Vo = 1 dE / w (l + z) (7) 
Aird 2 L dv 

and 

dR( z ) = X p ^l^-dz (8) 
1 + z dz 

In the above equations, di, = (1 + z)r is the distance luminosity, r is the proper 
distance, which depends on the adopted cosmology, dE gw /dv is the gravitational spec- 
tral energy emitted by a single source and v = (1 + z)v the frequency in the source 
frame. R c (z) is the cosmic star formation rate and X p is the mass fraction converted into 
neutron star progenitors with magnetic fields higher than 10 14 G. The (1+z) term in the 
denominator of eq. 8 corrects the star formation rate by the time dilatation due to the 
cosmic expansion. The element of comoving volume is given by 

with 

E(Qi,z) = {fl m (l + zf + n v } 1/2 (10) 

where tt m and Q v are respectively the present values of the density parameters due to 
matter (baryonic and non-baryonic) and the vacuum, corresponding to a non-zero cos- 
mological constant. A "flat" cosmological model (f2 TO + Cl v = 1) was assumed. Porciani 
& Madau (2001) provide three models for the cosmic star formation rate (SFR) history 
up to redshifts z ~ 5. Differences between these models are mainly due to the various 
corrections applied, in particular those due to extinction by the cosmic dust. In our 
computations, the three possibilities were considered. Since the final results are not sig- 
nificantly different according to the adopted SFR history, we present here only results 
derived from the second model, labelled SFR2. 

In our calculation, we have taken TO = 0.30 and Q v = 0.70, corresponding to the 
so-called concordant model derived from observations of distant type la supernovae 
(Perlmutter et al. 1999; Schmidt et al. 1998) and the power spectra of the cosmic mi- 
crowave background fluctuations (Spergel et al. 2003). The Hubble parameter H was 
taken to be 65 kms _1 Mpc _1 . 

Combining the equations above one obtains: 

o r \ SttG f Zsuv dE gw R SF R2{z) , .... 

or numerically: 

n gM = 7.5 x 10-s% 5 %v o £ UP ^^f^ 2 dz (12) 
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In the above equation, the SFR (in MQMpc _3 yr~ 1 ) as given by Porciani & 
Madau (2001) was derived for a matter dominated cosmology (f2 m = 1) with H = 
65kms _ Mpc - , but eqs. (11) and (12) already include the correction for the adopted 
cosmology. The upper limit of the integral is determined by the frequency cutoff in the 
source frame, namely 

Zsup = - 1 (13) 

Here we have assumed that the minimum rotation period is ~ 0.8 ms, which corresponds 
to a maximum gravitational wave frequency v sup w 2500 Hz. The cosmic SFR is quite 
uncertain for z > 5, but if one imposes the restriction z sup ~ 5, there are no practical 
consequences since objects formed at higher redshifts do not contribute significantly to 
the integral in eq. (12). 

3.2. Magnetar parameters 

The mass fraction of formed stars giving origin to magnetars is X p = x Jio° £(m)dm, 

where £(m) oc to~ 2 - 35 is the initial mass function, supposed to be given by Salpeter's 

law, and NS progenitors are supposed to form in the mass range 10-40 M . \ 1S the 

fraction of NS born with magnetic fields B > 10 14 G, which from our simulations is 

about 8%. Thus, X p = 3.9 x lO^M" 1 . 

The total spectral gravitational energy emitted by a rotating neutron star is given by 

dE gw _ dE gw | dt | 
dv dt dv 

where the frequency evolution, according to the analysis of Section 2.2, is controlled by 
the magnetic breaking. In this case, one obtains (RPOlb) 

dEgw 3 512(2 6 2 2 r 3 , . 

-&r = Kv ~ -^r* E h *if ■ (15) 

Notice that here we have corrected the above equation by the factor of two missing in 
the original equation derived by RPOlb. The ratio between the initial period Po and the 
magnetic braking time-scale tq is 

MS- 

For magnetars, the ellipticity is given by eq.(l) and in this case the constant K is given 
by 

K=^9 2 IL(<^/Bh i n*(a) (17) 
or, numerically 

K = 2.0x W 31 g 2 ll 5 R 2 10 M^B 2 14 sm 4 (a) (18) 
where I zz is given in units of 10 45 gem 2 . 
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Fig. 5. Density parameter £l gw for the stochastic background produced by magnetars for 
the two models considered in this work. Model A, solid curve and model B, dotted curve. 

According to our simulations, the average value of the quantity < (Bf 4 sin 4 a) > for 
the galactic population is 32.2, which will be assumed to be also representative of extra- 
galactic magnetars. Notice that if the migration of the magnetic dipole is introduced as 
in RP01, the average < (_Bi 4 sin 4 a) > should be replaced by < (i?^ 4 sin 2 a) > in eq. 17. 
However, this alternative model does not change significantly the numerical value of K, 
since the migration scale time is rather short (< 10 4 yr) and, as a consequence, most of 
the objects will be quickly found in an orthogonal configuration. 

4. Results and detectability 

Fig. 5 shows the density parameter Q gw as a function of the observed frequency for 
two NS models. Model A considers a normal interior (non-superconducting), having a 
representative deformation parameter g = 13 (Konno, Obata & Kojima 2000); model B 
assumes the existence of a type I superconducting interior, characterized by a deformation 
parameter g — 520 (Bonazzola & Gourgoulhon 1996). 

The density parameter Q gw reaches a maximum around 1200 Hz and has a high 
frequency cut-off at <~ 2500 Hz. The maximum value of the density parameter is about 
6 x 10~ 13 in the case of model A and ~ 10 -9 if the NS has a superconducting interior. 

If f is the duration of the signal produced by a single source, the stochastic back- 
ground from magnetars is expected to be continuous since the duty cycle satisfies 
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Table 1. Expected Q m in for the main detector pairs, corresponding to a flat background 
spectrum and an integration time T = 1 year, a detection rate a — 90% and a false 
alarm rate 7 = 10%. LHO and LLO stand for LIGO Hanford Observatory and LIGO 
Livingston Observatory 





LHO-LHO 


LHO-LLO 


LLO-VIRGO VIRGO-GEO 


initial 


4 x irr 7 


4 x 1(T 6 


8 x 1(T 6 8 x 1(T 6 


advanced 


6 x irr 9 


1 x 1(T 9 





D = dR(z)f (1 + z) » 1. Because it obeys Gaussian statistics and can be con- 
founded with the instrumental noise background of a single detector, the optimal detec- 
tion strategy is to cross correlate the output of two (or more) detectors, assumed to have 
independent spectral noises. The cross correlation product is given by (Allen & Romano, 
1999): 



Y= / sY(f)Q(f)s 2 (f)df (19) 

J —OO 

where 

is a filter that maximizes the signal to noise ratio (S/R). In the above equation, P\(f) 
and i-2(/) are the power spectral noise densities of the two detectors and T is the non- 
normalized overlap reduction function, characterizing the loss of sensitivity due to the 
separation and the relative orientation of the detectors. The optimized S/N ratio for an 
integration time T is given by (Allen, 1997): 



l )2 = 9_Hl r rxmiuf) 

'N> 8^ J J fe Pl (f)P 2 (f) 



& = ^TT df ^'Z ' (21) 



In the literature, the sensitivity of detector pairs is usually given in terms of the 
minimum detectable amplitude for a flat spectrum (Q gw equal to constant) (Allen & 
Romano, 1999) 

_ , — j. -i/ \ ^-x,^r M r2 (/) 



,<•*-'(»■> --*-<»r»ij[ «j4mm r ''' 2 <22) 



The expected minimum detectable amplitudes for the main detector pairs, after one 
year integration, are given in Table 1, for a detection rate a = 90% and a false alarm 
rate 7 = 10%. 

The power spectral densities used for the present calculation can be fond in Damour 
et al. (2001). 

O min is of the order of 10~ 6 — 10~ 5 for the first generation of interferometers combined 
as LIGO/LIGO and LIGO/VIRGO. Their advanced counterparts will permit an increase 
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Fig. 6. Integrand of the signal to noise ratio (see eq. 21) for advanced LIGO pairs: LHO- 
LHO (continuous line) and LHO-LLO (dot line) 



pair formed by the co-located and co-aligned LIGO Hanford detectors, for which the 
overlap reduction function is equal to one, is potentially one order of magnitude more 
sensitive than the Hanford/Livingston pair, if instrumental and environmental noises 
can be removed. However, the spectrum of magnetars is not flat and the maximum 
occurs out of the optimal frequency band of ground based interferometers, which is 
typically around few hundreds Hz, as in Fig. 6. Considering the co-located and co-aligned 
advanced LIGO interferometer pair, we find a signal-to-noise ratio S/N ~3x 1CP 6 and 
~ 5x 10~ 3 for models A and B respectively. A possible solution to increase the sensitivity 
at higher frequencies would be to use the interferometers in a narrowband (rather than 
in a wideband) configuration, as has been proposed for GEO 600, with the frequency 
of maximal sensitivity centred around 1 kHz. This could be achieved by choosing a 
high (rather than low) reflectivity of the signal in the recycling mirror, which would 
concentrate the power in a narrow band rather than distribute it in a wideband. The 
resonance frequency of the recycling cavity is tunable to the desired frequency by shifting 
the signal-recycling mirror. 

Unless the deformation parameter g is substantially higher than the present expec- 
tations, our results indicate that its contribution to the gravitational background is not 
detectable by the present generation antennas and also probably by their advanced coun- 



of two or even three orders of magnitude in sensitivity (CI. 
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Fig. 7. Theoretical limits for some cosmological gravitational background models and 
present estimates for the magnetar contribution. Other astrophysical contributions (in 
grey) are plotted for comparison: continuous background (continuous line) such as NS/NS 
coalescence (Regimbau & de Freitas Pacheco 2005) or rotating pulsars (updated from 
RPOlb) and shot noise (dot line) such as BH ringdown and bar mode due to the 
Maclaurcn/Dedekind transition in young and fast rotating NS (updated from Regimbau 
2001). The expected sensitivity resulting from the combination of two advanced interfer- 
ometers is also shown for comparison. 

terparts. A detection could be possible if the interior of the NS is constituted by a type 
II superconductor. Effective quark-quark interactions in the deconfined core may form 
a diquark condensate, which behaves as a superconductor of the second kind Blaschke, 
Sedrakian & Shahabasyan 1999), able to develop quite high fields and to produce a con- 
siderable magnetic distortion (g > 10 3 ). The deformation of the star can be increased if 
the magnetic field has a geometry other than dipolar, since higher order multipoles will 
also contribute to the mass quadrupolc. This is also true if the toroidal component is 
considered since its strength may be a factor of 10 higher than the poloidal component. 
Clearly these last possibilities require further work on NS structure when magnetic fields 
with complex geometries are considered. 

The background from magnetars may not be detectable as a signal by the first gen- 
erations of interferometers but it may behave as noise in the cosmological background 
frequency domain. 
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In Figure 7, we have plotted some theoretical predictions for the CGB (Allen 1997; 
Maggiore 1999; Buonanno 2004), our estimates for the magnetar contribution along with 
other AGBs and the expected sensitivity of two correlated advanced interferometers. The 
amplification of vacuum fluctuations at the transition between the radiation and matter 
eras, first discussed by Grishchuk (1974 & 1993) and Starobinsky (1979), is expected 
to produce a GW background whose spectrum and amplitude depend strongly on the 
fluctuation power spectrum developed during the early inflationary period. 

The standard de Sitter inflation predicts a spectrum that decreases as l// 2 in the 
range 3 x 10~ 18 — 1 x 10~ 16 Hz and then remains constant in a very large band up 
to frequencies of the order of MHz. The Cobe experiment provides an upper bound of 
Qgw = 9 x 10~ 14 (Maggiore 1999; Buonanno 2004) for the flat region, which is below 
the estimated background from magnetars in the frequency range 400 — 2500 Hz for 
model A and 70 — 2500 Hz for model B. These predictions have been questioned recently 
by Grishchuk (2005). In the more realistic scenario of "slow roll down" inflation, the 
predicted amplitude is even lower, by almost two orders of magnitude. A more interesting 
case arises from the string cosmology since the expected amplitude is around £l gw — 10~ 7 , 
well above the expected contribution from magnetars. However, other choices of the model 
parameters may lead to substantially lower amplitudes. Cosmic strings are also able to 
produce a CGB in the frequency range of the magnetar contribution, since the vibrations 
of these topological defects are sources of GW (Buonanno 2004). The spectrum is almost 
flat in a large frequency band (10~ 8 — 10~ 10 ) with an upper limit for the amplitude 
of 10~ 9 — 10~ 8 according to numerical simulations of a cosmic strings network that 
satisfy the bound imposed by the timing of pulsars. As it can be seen from fig. 7, the 
magnetar contribution from model B is barely comparable to the predictions of this 
model at frequencies around 1 kHz. However, according to a recent study by Damour 
and Vilcnkin (2000, 2001, 2005) the background from a network of cosmic strings at the 
end of Brane inflation may be not continuous but rather a collection of bursts emitted 
from cusps of oscillating string loops. In this case, it would be distinct in its statistical 
properties from the magnetar contribution. 

5. Conclusions 

In the present work, the contribution of magnetars to the stochastic background of gravi- 
tational waves was investigated, under the assumption that these highly magnetized stars 
are typical NS born in the high side tail of the initial field distribution (RPOla). Their 
geometric deformation is a consequence of such fields. 

The statistical properties of the magnetar population such as the magnetic field and 
cllipticity distributions and their relative orientation were derived by simulations and 
then by supposing that the extra-galactic population has similar properties. Our new 
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population synthesis, including new pulsars discovered by the Parkes multibcam survey, 
indicates that magnetars (neutron stars born with fields higher than 10 G) constitute 
8% of the total population. 

Two models for the neutron star interior have been considered in our calculations: 
in the first (model A), the NS has a normal conducting interior, leading to deformation 
parameters g in the range 1-20. Here, the value g = 13 was adopted (Konno, Obata & 
Kojima 2000). The second possibility (model B) considers that the NS develops a type 
I superconducting interior with a typical deformation parameter g — 520 (Bonazzola & 
Gourgoulhon 1996). 

We found that magnetars produce a continuous background (duty cycle >> 1) with 
a maximum around 1.2 kHz. This value is slightly dependent on the assumed maximum 
rotation velocity of the NS at birth, here taken equal to 1250 Hz (corresponding to a 
minimum initial rotation period of 0.8 ms). The different cosmic star formation rates do 
not considerably affect the resulting density parameter Q gw , whose amplitude depends 
on the deformation parameter <?, e.g., on the conducting properties of the NS interior. 
Amplitudes higher than the maximum found in this work, e.g., Ct gw ~ 10~ 9 , could be 
expected if neutron stars have deconfined cores and a diquark condensate forms a second 
kind superconducting medium. However, estimates of the resulting deformation parame- 
ter under these conditions remain to be done. In spite of the maximum amplitude being 
inside the frequency band of ground based interferometers, the expected signal is well 
below the sensitivity of the first generation of detectors, but exceeds the upper bound 
derived from COBE data. Comparing our estimates to the most optimistic expecta- 
tions of the gravitational background generated by some possible models, we find that 
the magnetar contribution does not overwhelm the relic emission expected from string 
cosmologies, but could be comparable around 1.0 kHz to the limits estimated for the 
emission from oscillating cosmic strings. Above 70 — 120 Hz the magnetar contribution 
surpass the theoretical limits of the inflationary CGB, suggesting that searches should 
remain well below those frequencies. In spite of the maximum amplitude lying inside the 
frequency band of ground based interferometers, the expected signal is well below the 
sensitivity of the first generation of detectors. 
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